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report 

-  S474V  is  «tied  by  the  <1  eophys  ic  s  Re  search  ^ 

Directorate  of  the  AltJf^OTce  Cambridge  Researoh  Gent^^r'^ob* 
^^taia>dffita  rnncegfriogothe  radiant  intensity  of  the  infrared  sky 
^  background  and  tho  amplitude  of  the  sky  gradient.  Measurement 
of  the  Infrared  sky  background  has  been  a  more  or  less  continuous 
program  since  July  1953  and  a  continuous  record  of  the  work 
accomplished  can  be  gleaned  from  as^samssk:^^  the  reports  issue^ 
which  are  included  in  the  reference  list.y  The^4a-ggeei  asnount-ef- 
m  the^o^  program  has  bee«i-4aken^nring  the 


.pxfiaant-coiitract;  AFir9^  604^-^73^;  Both  the  xadiometric  data  and 
the  gradient  data  discussed. ,  The  actual  data  have  been 

presented  in  B.eports  No.  1,  2,  and  3  preceding  this  final  report 
and  the  discussions  based  on  those  data.^  An  empirical 

form  has  been  worked  out  for  the  brightness  (fittt  taken  at  Key 
West,  i  lovida,  in  June  1959.  It  had  been  h(^ed  that  the  measure¬ 
ments  at  Key  West  might  result  in  a  mathematical  model  for  sky 
brightness  in  the  infrared.  This  might  ha^e  been  possible  had 
completely  cloudless  blue  skies  presen^d  themselves  for  meas¬ 
urement.  This  condition  never  occur j4d;  however,  work  has  been 
done  with  the  data  taken  to  result  in^  fairly  simple  and  possibly 
quite  u»?eful  equation  that  may  pro^,  with  further  investigation 

and  consequent  modification,  to  l/e  more  general  than  a  mathe- 

/ 

matical  model  applying  only  to^lue  sky. 

Gradient  data  is  also  meussed  and  some  work  has  been 
done  to  determine  the  valuM  of  the  most  likely  gradients  for  blue 
sky,  clouds,  and  cloud  ewes.  A  theoretical  mathematical  model 
for  the  infrared  backgr^nd  gradient  is  presented  which  includes 
some  considerations  TOncerning  the  spectrum  expected  viewing 
the  earth  Lackground^from  various  altitudes.  Some  discussion 
is  included  concerning  the  use  of  gradient  data  in  the  design  of 
equipment.  ^ 


This  final  report  also  includes  a  sumn-  Jiry  of  the  background 
investigwt  on  work  from  its  inception  in  193i  with  the  building  of 
the  instrumentation  used  to  obtain  the  data*  Some  of  the  more  im¬ 
portant  phases  of  the  work  are  discussed  in  some  detail  so  that 
this  report  will  bring  most  of  the  past  work  into  one  bound  volume. 
The  theoretical  mathematical  model  for  the  background  gradient 
gives  indications  of  areas  for  future  measurement  work  and  this 
is  discussed  in  the  summary. 

Contract  AF  19(60^1) -347  3  has  been  performed  by  the  Infrared 
Dep  rtment  of  Ramo- Wooldridge «  a  Division  of  Thompson  Ramo 
Wooldridge,  Inc.  Personnel  directly  connected  with  the  perform¬ 
ance  of  this  Contract  are: 


T.  R.  Whitney 
G.  F.  Aroyan 
P.W.Shadle 
.R.  £.  Eisele 
J.W.Fisk 
C«  J*  Somol 


Member  Senior  Staff 
Member  Senior  Staff 
Member  Technical  Staff 
Member  Technical  Staff 
Member  Technical  Staff 
Research  Assistant 


Department  Head 
Asst.  Dept.  Head 
Group  Head 
Project  Head 


DISCUSSION  OF  KEY  WEST  BACKGROUND  DATA 


Tlie  data  published  in  Reports  No.  Z  and  3  were  taker  at 
Key  West,  Florida,  in  June  1959*  as  part  of  a  measurement  field 
trip  in  cooperation  IRMP  59/60.  The  data  was  obtained  with 
a  combined  filter  radi  meter  and  gxa  lient  meter,  the  design  and 
operation  of  which  is  t  mpletely  described  in  previous  reports. 

The  calibrations  and  corrections  used  in  reduction  of  the  data  are 
included  in  the  two  reports  (4)  containing  the  ts.bulations.  The 
^urpose  in  this  section  will  be  to  summarize  the  data. 

Key  West  is  at  24^  34'  latitude  and  81^  45'  longitude.  It  is 
thus  just  north  of  the  Tropic  of  Capricorn  and  well  within  the  hur¬ 
ricane  belt.  On  J\me  21,  the  sun  is  at  the  highest  elevation  and 
passes  very  close  to  zenith  at  noon  solar  time  and  also  follows  a 
general  East- West  azimuth.  The  plot  of  sun  ordinates  with  re¬ 
spect  to  local  standard  time.  Figure  1,  shows  that  the  path  of  the 
sun  does  not  vary  more  than  5  degrees  from  this  azimuth  between 
the  hours  of  10  a.  m.  and  3  p.  m.  and  no  more  than  20  degrees 
through’ithe  course  of  the  day.  It  was  hoped  that  radiometric  or 
brightness  data  taken  on  cloudless  skies  would,  for  this  reason, 
lend  itself  more  easily  to  evaluation  and  lead  to  a  simple  function 
for  sky  brightness  in  the  infrared,  thus  making  conclusions  more 
general  for  the  area.  Simply  because  of  the  location  and  the  time 
of  year,  the  mixing  of  the  air  mass  becomes  quite  turbulent  and 
leads  to  constant  cloudiness,  recurrent  storms  and  brisk  winds. 
The  data  is  indicative  of  these  conditions  and  represents  the  ra¬ 
diative  character  of  the  sky  only  for  this  location  and  time  of  year. 

RADIOMETER  DATA 

Data  taken  June  16th  and  June  17th  arc  for  completely  over¬ 
cast  skies  which,  to  the  eye,  were  completely  gray  and  feature¬ 
less.  The  sun  position  was  indicated  by  a  large  bright  area  in  the 
the  cloud  cover  due  to  heavy  scattering.  Total  radiation  under 


such  conditions  was  measured  using  the  bolometer  with  no  spec¬ 
tral  filters  in  the  system.  The  bolometer  spectral  range  is  limited 
by  a  sodium  chloride  window  and  represents  a  wavelength  range  from 
0  to  16  microns.  Polar  plots,  h  igure  2  and  3,  sum  the  data  points. 
The  sun  position  is  indicated  by  an  asterisk  and  represents  the  mid¬ 
point  of  the  sun's  travel  through  the  course  of  the  measurements. 
About  30  minutes  is  required  to  gather  the  information  for  one 
plot  during  which  time  the  sun  would  have  moved  a  total  of  7.  5 
degrees.  The  two  plots  are  quite  similar  with  small  local  differ¬ 
ences  which  are  no  doubt  due  to  differences  in  thickness  of  cloud 
cover.  The  energy  levels  in  Figure  2  are  slightly  lower  near  the 
sun  than  for  Figure  3,  and  for  an  earlier  time  of  day,  but  the  dif¬ 
ferences  are  not  great  and  are  no  doubt  due  to  changing  conditions 
in  overcast.  A  low  area  is  indicated  in  Figure  3  about  60  degrees 
west  of  the  sun  position,  but,  since  there  is  no  indication  of  this 
in  the  previous  day's  run,  the  cause  again  is  probably  differing 
thickness  in  cloud  cover.  Since  this  spectral  region  includes  the 
8  to  12  micron  "window”,  these  differences  may  be  due  to  warmer 
and  colder  cloud  patches  with  the  cooler  areas  being  thicker.  This 
conjecture  is  more  plausible  under  such  conditions  than  to  presume 
the  cause  to  be  solely  scattering  in  the  short  wavelength  region 
since  the  sky  was  visibly  uniformly  gray  with  consequent  low  illu¬ 
mination  levels  as  borne  out  by  filter  measurements  in  the  near 
infrared  under  the  same  conditions. 


Figure  4  is  a  detailed  plot  of  the  overcast  sky  at  4.  25  mi¬ 
crons.  This  information  was  obtained  immediately  following  the 
total  radiation  measurement  in  Figure  3,  but  no  particular  corre¬ 
lation  is  indicated  between  the  two  plots.  A  much  more  uniform 
ofet  of  curves  probably  would  have  resulted  had  the  spectral  filter 
been  narrower.  The  filter  at  4.  25  microns  has  a  one  micron 


bandwidt^  and  consequently  radiation  is  admitted  on  both  sides  of 
the  CO^  absorption  band.  Several  high  energy  regions  are  indi¬ 
cated  on  the  diagram  approximately  75  degrees  from  the  sun  and 


an  increasing  gradient  is  indicated  in  the  southeast  area  down  to 


the  horizon.  The  measurements  in  this  direction  were  over  water 
which  extended  perhaps  one -half  mile  and  directly  beyond  this  were 
concr^^^  runways  towards  the  end  of  the  airstrip  at  Boca  Chica 
(Naval  Air  Station  at  Key  West).  Probably  because  of  the  inclem¬ 
ent  weather,  there  were  no  aircraft  operating  on  either  of  these 
two  days.  The  total  radiation  map  does  in  a  general  way  show  a 
high  tendency  in  this  same  area.  Several  days  later,  June  23rd, 
it  was  possible  to  obtain  data  on  blue  sky  patches  in  a  very  cloudy 
sky.  Figure  5  presents  this  group  of  measurements  for  compari¬ 
son  with  Figure  4.  In  general,  the  data  are  rather  similar.  The 
t}  "c  sky  observations  do  not  show  the  details  of  the  overcast  sky 
and  the  gradient  i^  much  more  uniform.  A  general  minimum  area 
can  be  seen  at  about  70  degrees  from  the  sun.  Jime  23rd  was  a 
bright  sunny  day  with  many  large,  white,  cumului  clouds.  The 
ambient  temperature  was  84^  F  on  June  17th  and  85"^  F  on 
June  23rd,  but  the  sun  was  at  a  higher  elevation  for  the  measure¬ 
ments  of  the  23rd.  The  values  at  the  horizon  are  higher  than  the 
blackbody  radiation  at  4.  25  microns  for  the  ambient  35^  F  which 
is  140  pw/cm^  *  steradian  •  micron.  Again,  this  may  be  due  to 
radiation  outside  the  CO^  absorption  band. 

Polar  plots  could  be  presented  for  the  other  spectral  regions 
except  that  the  details  are  so  numerous  that  the  plots  are  very  dif¬ 
ficult  to  prepare  and,  as  a  consequence,  would  be  difficult  to  in¬ 
ter^  ret.  It  is  felt  that  a  spectral  plot  showing  highs,  lows  and 
averages  for  the  day  give  a  much  better  general  idea  of  the  radia¬ 
tion  characteristics.  This  information  is  given  for  June  16th  in 
Figure  6  for  heavy  overcast.  The  average  readings  are  shown 
bracketed  with  a  shaded  area  designated  as  "measurement  accu¬ 
racy".  The  measurement  accuracy  is  taken  as  t  6%  since  both 
the  energy  calibration  and  the  voltage  calibration  each  repeat 
within  2%  and  the  data  is  read  from  the  record  probably  within 
/-%.  Whenever  the  low  values  were  read,  a  variation  of  about 
1 5|jLw/cm  •  steradian  •  micron  was  noticed  over  a  10  to  15 
second  time  duration.  This  has  been  attributed  to  variations  in 


the  background  since  it  does  not  occur  in  the  laboratory  when 
looking  into  the  blackbody,  nor  is  it  periodic.  Further,  no  cor-* 
relation  between  filtex  regions  could  be  observed,  i.e.  the  energy 
level  might  be  increasing  slightly  in  one  filter  region  while  de¬ 
creasing  in  another.  Motion  of  the  clouds  may  be  the  largest  con¬ 
tributing  factor  although  variations  of  this  small  amplitude  could 
be  due  to  many  factors  in  the  background.  The  variation  has  been 
*'iided  to  the  accuracy  envelop.  Average  spectral  data  for  June  17, 
i,  1  and  23  are  given  iu  Figures  7,  S  and  9.  The  measurements 
taken  June  19th  are  designated  as  '*cloudy-bright*'  which  means 
that  the  sky  was  overcast  but  clearing  with  a  consequent  increase 
in  the  brightness.  The  23rd  was  quite  cloudy  with  patches  of  blue 
sky  showing  all  over  which  gave  the  sky  a  mottled  appearance. 

All  data  for  this  day  are  for  blue  sky  under  this  condition. 

gradient  data 

The  gradient  data  has  been  sorted  into  three  general  classi¬ 
fications:  (1)  Blue  Sky,  (2)  Clouds,,  and  (3)  Clouds  and  Sky. 
Those  for  blue  sky  are  for  blue  patches  of  sky  under  very  bright 
cloudy  sky  conditions.  It  is  indicative  of  the  conditions  at  Key 
West  in  June  that  only  90  such  measurements  could  be  made  dur¬ 
ing  the  whole  measurement  period.  Since  a  6  degree  scanning 
circle  was  used  to  develop  the  background  function,  patches  of 
sky  smaller  than  10  or  12  degrees  were  ignored  because  of  the 
possibility  of  inadvertently  including  clouds  in  the  field.  To  en¬ 
hance  the  usefulness  of  the  gradient  data,  each  classification  was 
sorted  to  determine  the  probable  amplitude  occurrence  at  particu¬ 
lar  frequencies.  Histograms  were  constructed  for  ten- space  fre¬ 
quencies  from  3  to  150  waves  per  radian.  From  these,  probabil¬ 
ity  curves  were  developed  and  95%  distribution  limits  were 
determined.  As  an  example.  Figure  10  shows  how  the  measured 
amplitudes  for  3  waves  per  radian  were  distributed  for  the  blue 
sky  data.,  Figure  1 1  is  a  photograph  showing  the  distributions 


for  cloudft  erected  on  a  frequency -amplitude  chart.  This  treats 
ment  of  gradient  data  provides  95%  distribution  limit  charts  for 
the  three  general  classifications.  Figures  12,  13,  and  14  repre¬ 
sent  the  complete  statistical  summary  of  the  sky  gradient  in  the 
lead  sulfide  region  for  Key  West  during  the  month  of  June.  These 
charts  include  all  positional  points  measured  and  for  this  reason 
are  completely  general.  These  charts  are  presented  with  the 
energy  ordinate  squared  and  are^  therefore,  the  one  dimensional 
power  spectra  for  the  three  sky  conditions.  The  center  solid  line 
is  the  peak  of  the  probability  curves  with  the  upper  and  lower 
solid  lines  determining  the  upper  and  lower  limits  between  which 
95%  of  the  total  number  of  measurements  lie.  Dashed  lines  of 
various  slopes  are  included  for  comparison.  Similar  measure¬ 
ments  made  from  Mt.  Wilson  in  1956  showed  similar  character¬ 
istics  with  the  power  spectra  running  slopes  to  -3. 

Because  secondary  peaks  regularly  appear  in  the  probabil¬ 
ity  curves  for  clouds,  these  are  presented  in  Figures  15,  16,  and 
17.  A  definite  trend  to  a  secondary  slope  at  higher  amplitude 
levels  is  indicated  for  the  cloud  measurements.  No  explanation 
of  this  is  attempted  here  except  to  point  out  that  the  classification 
"Clouds*'  is  quite  general  and  further  breakdown  of  this  classifi¬ 
cation  into  cloud  types  might  reveal  more  detailed  trends. 

At  times,  gradients  were  measured  on  the  blue  sky  to  look 
for  trends  from  horizon  to  zenith  for  one  sun  position.  No  sig¬ 
nificant  trend  was  observed.  Figures  18  and  19  can  be  considered 
as  representative  of  this  type  of  data.  Figure  18  shows  the  spectra 
observed  for  different  elevations  for  the  azimuth  plane  nearly  at 
right  angles  to  the  plane  of  the  sun,  while  Figure  19  gives  similar 
information  for  the  azimuth  plane  of  the  sun. 
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Figure  2.  Overcast  Sky,  June  16,  1959.  Key  West,  0-16  Microns. 
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Figure  3.'.  Overcast  Sky,  Jiine  17,  1959 •  Key  West,  Microns. 
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Figure  5.  Blue  Sky  Areas  in  a  Cloudy  Sky,  June  23,  1959»  Key  West, 
4.  25  Microns. 


Figures.  Average  Spectral  Distribution,  Overcast  Sky, ‘June  16,  1959. 
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Figure  8.  AverJige  Spectral  Dietribution,  Overcast  Sky,  Xloudy -Bright", 
June  19i  1959* 

'  15 


Figure  9*  Average  Spectral  Distribution,  Blue  Sky  Areas  in  a  Cloudy  Sky 
June  23,  1959. 
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Figure  10.  Amplitude  Probability  for  Blue  Sky  at  3  Waves  per  Radian. 
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Figure  12!.  Ninety -Five  Percent  Distr  ibution  Limits  for  Blue  Sky  One 
Dimensional  Power  Spectrum. 
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Figure  13.  Ninety -Five  Percent  Distribution  Limits  for  Clouds  One 
Dimensional  Power  Spectrum. 
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Figure  14.  Ninety -Five  Percent  Distribution  Limits  for  Cloud  and  Sky 
One  Dimensional  Power  Spectrum. 
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Figure  15.  Individual  Amplitude  Probability  Curves -Clouds. 
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Figure  19*  Individual  Spectra  for  Blue  Sky  for  the  Sun  Azimuth  Plane. 
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INITIAL  CONSIDERATIONS  UNDERLYING 
A  MODEL 

One  of  the  stated  objectives  of  the  Key  West  measurement 
field  trip  was  to  obtain  data  which  might  produce  a  mathematica.1 
model  for  sky  brightness.  What  might  be  classed  as  a ‘rigorous 
mathematical  model  did  not  culminate  from  these  data.  Sky  con¬ 
ditions  were  widely  variant  and  cloudless  sky  did  not  exist  during 
the  whole  period.  Nevertneless ,  a  great  amount  of  effort  has  been 
applied  to  make  a  beginning  attempt  and  future  effort  possibly  will 
result  in  a  usable  practical  model  for  the  infrared  sky  background. 
This  section  of  the  report  discusses,  in  general  fashion,  the  basic 
considerations  underlying  the  description  of  a  model.  The  next 
section  describes  some  effort  to  fit  the  sky  brightness  values  with 
simple  expressions  empirically. 

Figure  20  depicts  a  simple  model  of  the  earth  and  its  atmos¬ 
phere.  The  atmosphere  consists  of  a  mixture  of  permanent  gases 
that  exist  at  a  fixed  percentage  of  the  total  and  other  gases  whose 
presence  is  highly  variable.  The  accompanying  table  summarizes 
the  composition  of  the  atmosphere,  neglecting  those  gases  which 
are  present  in  only  small  aixvounts. 


Table  I 


Constituent 

Per  Cent 

p 

8, 

Nitrogen 

78 

^  ik 

m  a 

Oxygen 

21 

ri  S 

.r\  '  e 

^Carbon  Dioxide 

0.033 

8  E 

n 

t  y 

Argon 

0.  93 

o 

a 

'  i 

Water  ‘Vapor 

0-2 

1 

8 

a 

Ozone 

0-,  07  p.  p.  m. 

b  ® 

m 

(ground  level) 

1 

^5  ^ 

1-3  p.  p.m.  (20-30  KM) 
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Throughout  the  troposphere  and  the  stratosphere  the  perma¬ 
nent  gases  exist  at  nearly  the  same  percentage,  whereas  the  water 
vapor  content  shows  a  rapid  falloff  with  increasing  height.  Although 
water  vapor  is  generally  less  than  two  percent  of  the  atmosphere, 

It  absorbs  about  six  times  as  much  solar  radiant  energy  as  all  the 
other  gases,  and  also  accounts  for  almost  all  the  absorption  of  the 
earth's  radiation. 

In  addition  there  are  a  variety  of  particles  suspended  in  the 
atmosphere,  which  comprise  the  atmospheric  aerosol.  Most  im¬ 
portant  are  the  water  droplets  and  ice  crystals  of  clouds.  Other 
particles  include  dust,  smoke,  sea  salt  particles,  and  minute  or¬ 
ganisms  such  ds  pollen,  spores,  bacteria,  etc. 

The  next  table  indicates  the  depletion  and  absorption  of  inci¬ 
dent  solar  energy  under  clear  skies  and  cloudy  skies,  as  well  as 
the  mean  conditions.  (The  average  cloud  cover  of  the  earth  is  52 
percent. ) 

Table  II 


Clear  Skies 

Per  Cent 

Atmospheric  Scattering 

9% 

Surface  Reflection^ 

6% 

Total 

15% 

Absorption  by  Water  Vapor 

13% 

Absorption  by  the  Other  Constituents 

3% 

Total 

16% 

Total  Depletion 

31% 

Total  Absorption 

69% 

Overcast  Skies 

Reflection  and  Backs catte ring 

55% 

Absorption 

10% 

Total  Depletion 

65% 

Total  Absorption 

35% 
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Table  II  (Cont'd. ) 


Mean  Conditions  Per  Cent 

Absorbed  through  clear  skies  33% 

(.  43  X  .  69) 

Absorbed  und^P‘^  cloud  cover  18% 

(.  52  .  35) 


Total  Incident  Energy  Absorbed  51% 


As  far  as  best  observations  indicate,  the  surface  of  the 
earth  is  not  undergoing  a  change  in  temperature.  Therefore,  there 
must  be  a  balance  between  incident  energy  and  outgoing  energy. 

In  the  unoccluded  portions  of  the  earth,  radiation  is  transmitted 
through  atmospheric  "windows"  and  -  in  the  bands  of  water  vapor 
and  carbon  dioxide  absorption  -  energy  is  radiated  to  space  as 
blackbody  radiation  at  the  temperature  of  the  stratosphere.  When 
the  sky  is  overcast  the  only  emission  is  at  the  blackbody  tempera¬ 
ture  of  the  stratosphere.  It  is  known  that  these  effects  vary  over 
the  surface  of  the  earth,  with  season,  latitude,  region,  time,  etc. 
As  a  simplification  of  an  extremely  complex  process,  the  earth 
will  bv  assumed  to  be  uniform,  i.  e.  ,  each  area  will  be  assumed 
identical  with  equal  probability  of  being  occluded.  The  next  table 
presents  a  summary  of  the  earth -atmosphere -solar  radiation  heat 
balance.  It  is  assumed  that  100  units  of  radiation  are  incident 
upon  the  atmosphere. 

SOLAR  (Space)  BALANCE 

IN  OUT 

100  units  Solar  Radiation  scattered  upward 

Reflected  Solar  Radiation 

Emitted  Radiation  from  Absorbing 
Transparent  Region 

Emission  from  the  Earth  in  the 
Transparent  Region 


9 

27 

47 

17 

100 
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EARTH  SURFACE  BALANCE 


IN  OUT 


Direct  Solar  Radiation 

32 

Emission  from  the  earth  in 
atmospheric  absorption 

120 

regions. 

Solar  Radiation  scattered 

19 

Heat  carried  to  the  atmos  - 

23 

downward. 

106 

phere  in  the  hydrological 
cycle. 

Infrared  emission  from  the 

Emission  from  the  earth  in 

17 

atmosphere. 

Heat  transported  down  by 

3 

transparent  regions. 

eddy  currents. 

160 

160 

ATMOSPHERIC  BALANCE 

IN 

OUT 

Earth  emission  in  absorbing 

120 

Infrared  emission  to  space. 

47 

regions 

Heat  carried  to  the  atmos¬ 

23 

Infrared  emission  to  the 

106 

phere  in  the  hydrological 

earth. 

cycle. 

Solar  radiation  absorbed. 

13 

Heat  transported  downward 
by  eddy  currents. 

3 

156 

156 

The  foregoing  relatione  delineate  the  gross  phenomena  of  the 
model.  Subsequent  paragraphs  will  discuss  in  detail  seme  of  the 
important  processes. 

ABSORPTION 

Figure  21  presents  typical  absorption  spectra  of  the  gases 
important  in  the  absorption  process  of  the  atmosphere.  (The  scales 
are  in  terms  of  transmission;  the  upper  limit  being  100  percent  trans¬ 
mission,  the  lower  limit  0  percent  transmission  or  total  absorption.  ) 
The  solar  spectrum  is  included  as  a  reference.  It  can  be  seen  that 
the  product  of  the  three  transmission  curves  for  the  gases  considered 
closely  approximates  the  solar  spectrum,  indicating  that  they  are 


the  only  gases  that  need  be  considered  in  a  simple  model.  The 
model  indicates  that  in  the  lower  atmosphere  only  water  vapor  and 
carbon  dioxide  need  be  considered,  and  that  ozone  is  the  first  gas 
to  effect  solar  radiation.  Since  the  concentration  of  ozone  i«#  very 
low  at  low  altitudes,  and  only  one  percent  of  the  atmosphere  lies 
above  the  ozone  band,  these  approximations  should  not  introduce 
significant  errors.  In  the  region  between  the  earth  and  the  ozone 
band,  the  role  of  carbon  dioxide  becomes  more  important  with  in¬ 
creasing  altitude.  Near  ground  level,  as  has  been  indicated,  water 
vapor  accounts  for  the  majority  of  the  absorption,  however,  its 
concentration  falls  off  rapidly  with  altitude,  whereas  the  partial 
pressure  of  carbon  dioxide  relative  to  the  total  pressure  is  essen¬ 
tially  constant. 

Howard,  Burch,  and  Williams^  ^  ^  discuss  the  problem  of 
absorption  by  water  vapor  and  carbon  dioxide  in  considerable  detail. 
After  a  series  of  measurements  of  the  absorption  of  these  gases  at 
various  pressures  and  with  various  path  lengths,  it  was  found  that 
two  types  of  empirical  relationships  were  sufficient  to  represent  the 
total  absorption  of  any  band  in  the  0.  9p  to  15p  region: 


where 


jA^dv  =  cw^^^  (P  +  pJ*' 

jAydv  =  C  +  Dlogjo  w+Klog,Q(P  +  p) 
jA^dv  =  the  total  absorption  of  a  band 


(weak  absorption  bands) 
(strong  absorption  bands) 


w  =  total  absorber  concentration  (atom -cm  for  CO^, 
pr.  cm  for  K2O) 

P  =  the  total  pressure 

p  =  the  partial  pressure  of  the  absorber. 

c,  k,  C,  D,  K,  =  constants  for  a  given  band. 

The  reference  contams  tabulations  of  the  constants  in  the  above 
relationships.  For  a  particular  absorption  band,  there  is  a  value  of 
total  absorptioi  below  which  the  ’’weak'*  relationship  holds,  and  above 
which  the  “strong’-  relationship  is  applicable. 
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There  have  been  two  simple  theoretical  models  created  to 
explain  the  observed  absorptions  by  water  vapor  and  carbon  dioxide; 
the  Elsas 8 er  Model  and  the  Goody  Model.  It  has  long  been  realized 
that  a  spectral  "line"  must  have  some  finite  width  because  of  broaden¬ 
ing.  For  H^O  and  CO^i  under  the  temperature  and  pressure  of  the 
lower  atmosphere,  the  finite  widths  of  the  spectral  lines  are  du^  to 
molecular  ccllisions.  For  a  X>orentz -shaped  line,  which  is  the 
shape  appropriate  to  the  model  atmosphere,  the  absorption  coeffi¬ 
cient,  kv  is: 


where 


IT  (v  -  v^)  +  a 

S  =  line  strength 
a  =  half  width  at  maximum 
V  =  frequency 
v^  =  center  frequency 


The  line  intensity  S,  is  essentially  independent  of  temperature 
and  a  is  a  function  of  temperature,  T,  and  pressure,  P. 


where  the  subscript  refers  to  the  values  at  standard  conditions. 

The  detailed  structure  of  a  CO^  band  is  a  fairly  regular, 
approximately  evt^nly  spaced,  sequence  of  absorption  lines,  where 
adjacent  lines  have  nearly  equal  intensity.  In  contrast,  an  H^O 
band  consists  of  many  irregularly  spaced  lines  of  very  unequal 
intensities. 

Elsasser  treated  the  case  of  an  idealized  band  consisting  of 
an  infinite  number  of  equally  intense,  equidistant  Lorentz- shaped 
lines.  For  this  case  the  absorption  coefficient  -  assuming  lines  at 
V  =  o,  -  2d,  1 4d,  ---  is 
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kv  = 


a 

(v  -  nd)  +  a 


Letting  < 


2ttv  j  ,  Zira 

a  =  —  and  b  - 

d  d 


S  sink  b 


d  cosh  b  -  cos  a 


If  the  Incident  energy  is  constant  throughout  the  spectrum,  the 
absorption  A,  is: 


2ir 


1 . 


da 


If  the  line  spacing,  d,  is  large  compared  to  the  half -width,  c*,  then: 


kv  =  k(a)  = 


S  b 


2d  sin 


.  2|a\ 

It) 


and  A  =  erf 


N 


irSaw 


where 


erf{x)  =  — 
IT 


The  average  absorption  is 


A 


where  i  =  ^ 


For  a  frequency  interval  which  contains  several  approximately 
equally  intense,  equidistant  lines,  where  there  is  no  absorption  at 
the  center  of  the  band  from  lines  outside  the  interval: 


The  absorption  bands  of  CO^  are  sufficiently  regular  to  satisfy 
the  requirements.  Thus  the  Elsasser  model  should  apply  to  CO^. 
Experiments  have  demonstrated  that  there  is  good  agreement  be¬ 
tween  measured  and  predicted  absorption. 
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Goody  developed  a  statistical  or  random  model  for  a  disordered 
band,  which  has  been  shown  to  fit  H^O  absorption  quite  well.  He 
considers  the  fractional  transmission  at  the  center  of  a  frequency 
interval  nS  wide  where  n  is  the  number  of  lines  of  Lorentz  shape 
in  the  interval,  and  6  the  mean  line  spacing.  The  line  spacings 
are  assumed  to  be  random,  and  the  line  intensities  which  are  un- 
correlated  with  position,  are  assumed  to  follow  the  distribution 
law  P(S) 

With  these  assumptions,  it  is  possible  to  show  that  the  trans- 

\ 

mission  at  the  center  frequency  of  a  thickness  w  of  absorbing  gas 
is: 

T(a,w^)  =  exp 


Howard,  Burch,  and  Williams  chose  to  represent  the  variation  of 
(7  over  a  band  by  the  relation  (7  =  Then,  letting  )3  =  k/6 

and  y  =  k/wa. 


Values  of  the  parameters  |3  and  y  were  determined  from  experi¬ 
mental  data,  and  a  final  result  obtained. 


where  w^  is  absorber  concentration  required  to  drop  the  trans¬ 
mission  to  50  percent  at  a  given  frequency.  This  relation  was 
found  to  fit  all  the  H2O  bands  considered  in  the  referenced  report. 
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SCATTERING 


In  addition  to  absorption  by  the  atmospheric  gases,  energy 

incident  upon  the  unoccluded  atmosphere  is  scattered  by  the  air 

molecules  and  by  the  atmospheric  aerosol. 

For  ’‘clear*'  air,  i.  e.  ,  assuming  that  there  are  no  clouds  or 

haze  particles  present,  the  Rayleigh  theory  of  scattering  by  air 

/  7) 

molecules  as  extended  by  Tousey  and  Hulburt'  '  can  be  used  to 
calculate  the  amount  of  both  primary  and  multiply  scattered  energy 
reaching  an  observer.  Because  the  referenced  report  contains  an 
extensive  development  of  the  theory,  it  will  only  be  summarized 
here.  Their  notation  is  as  follows: 

1.  the  zenith  angle  of  the  sun. 

2.  I,  the  zenith  angle  of  observation 

3.  (I,  the  scattering  angle. 

4.  Z,  projection  of  0  in  the  horizontal  plane. 

cos  0  =  cos  4  cos  1^4  ^  cos  Z 

5«  1  ,  intensity  of  sunlight  outside  the  atmosphere  as  a 

function  of  wavelength,  K 

6.  6i  the  polarization  defect  -  the  ratio  of  the  weak  to  the 

strong  polarized  component  of  sky  light  at  0  =  tt/Z 
(It  is  assumed  to  have  the  value  6  =  0.  04 

7.  r,  diffuse  earth  reflectivity  (abedo).  (An  average  val:- 

value  of  0.  20  is  assumed  for  model  earth) 

8.  index  of  refraction  for  air  -  function  of  wavelength. 

9.  n,  the  molecular  density  of  the  clear  atmosphere  (mole- 
^  cules  /cm^) 

10.  X,  thickness  of  the  atmosphere  above  the  point  of  obser- 

» 

vation,  reduced  to  NTP. 

Scattered  light  reaching  an  observer  consists  of  light  scattered 
directly  (primary  scattering)  and  light  which  is  scattered  more  than 
once  (mulitple  scattering).  The  total  observed  is  the  sum  of  the 
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Because  of  the  cos^  0  term,  the  scattering  should  be  maximum 
in  the  direction  of  the  incident  light  (0  =  0^)  and  in  the  opposite 
direction  (0  =  180^)  and  a  minimum  in  the  plane  at  right  angles  to 
the  incident  light  (0  =  90°).  The  polarization  of  the  light  should  be 
complete  in  this  plane  and  zero  at  0  =  0°»  0  =  180°.  In  practice  - 
because  of  multiple  scattering,  the  polarization  will  be  more  com¬ 
plicated.  The  attenuation  coefficient  for  pure  air  has  a  theoretical 
value  of  ^  =  0,  01 26 /km.  However,  when  theoretical  predictions 
were  compared  with  values  measured  above  10,000  feet,  it  was 
found  that  ^  had  to  be  increased  to  =  0.  017 /km  to  get  good  agree¬ 
ment.  Below  10,000  feet  scattering  by  dust,  larger  haze  particles, 
etc.  ,  was  so  large  as  to  make  the  assumption  of  a  Rayleigh  atmos¬ 
phere  unwarranted.  The  fact  that  had  to  be  increased  from  the 
pure  air  value  indicates  that  there  are  a  significant  number  of 
haze  particles  at  higher  altitudes.  However,  a  comparison  between 
measured  values  and  values  calculated  with  =  0.  017  /km  (presented 
in  the  reference),  shows  that  the  correction  is  sufficient  to  enable 
the  theory  to  give  useful  results. 

In  the  region  below  10,000  feet  the  most  significant  contribution 
to  scattering  is  due  to  the  particles  of  the  atmospheric  aerosol.  The 
scattering  of  incident  light  by  the  particles  of  the  aerosol  depends 
upon  their  size,  number,  index  of  refraction,  and  the  wavelength 
of  the  light.  The  most  important  class  of  particles  is  liquid  droplets. 
As  far  as  is  known,  the  only  possible  mechanism  for  the  creation 
of  liquid  droplets  from  water  vapor  depends  upon  the  presence  of 
condensation  nuclei.  The  decrease  in  visibility  when  the  relative 
humidity  exceeds  about  70  percent  indicates  that  these  nuclei  are 
hygroscopic.  Measurements  by  Junge'  '  indicate  that  the  distri¬ 
bution  of  nuclei  sizes  follows  this  law. 

N  =  K/R^ 

where  N  =  number  of  particles  per  unit  volume 
K  =  constant 
R  =  particle  radius. 


This  law  appears  to  be  valid  in  the  region  R  =  0.  05p  to  R  =  lOjio 
The  maximum  particle  concentration  lies  between  0.  01^  and  0.  05(i; 
for  smaller  R  the  distribution  decreases  rapidly  with  decreasing 
radius.  On  the  average,  the  law  applies  over  the  earth’s  surface. 

The  condensation  nuclei  grow  by  condensing  water  on  the  surface. 

As  the  droplet  grows  the  solution  concentration  decreases  and  the 
index  of  refraction  chainges.  Since  there  is  a  wide  range  of  possible 
particle  sizes,  particle  scattering  can  vary  widely. 

The  original  solution  of  the  problem  of  scattering  by  spherical 
particles  is  due  to  Mie.  Considering  scattering  alone,  the  expression 
for  the  transmission  of  light  follows  an  exponential  law: 

T=e-^ 

where  T  =  transmission 

A  =  scattering  coefficient 
X  =  path  length 

Using  the  techniques  of  classical  electromagnetic  theory,  Mie 
showed  that  the  scattering  coefficient  could  be  expressed  as: 

A  =  ir  NK 

where  r  -  particle  radius 

n  =  number  of  particles /cm^ 

K  =  the  scattering  area  ratio 

The  number  K  relates  the  effective  scattering  ratio  to  the 
geometrical  cross  sectional  area.  It  is  a  function  of  the  index  of 
refraction  of  the  particle  and  the  ratio  of  particle  radius  to  wave¬ 
length.  Figure  22  is  a  detailed  graph  of  the  variation  of  K  with 
the  parameter  a  =  2TTr/X  for  a  pure  water  droplet  (m  =  1.  33). 

Figure  23  presents  K  for  several  values  of  the  index  of  refraction. 

(m  1.  55  corresponds  to  dry  salt,  which  is  believed  to  be  the  most 
important  condensation  nucleus. )  The  curves  for  any  value  of  m 

tend  to  K  =  2  with  large  a-  For  low  values  of  CLt  corresponding  to 

-4 

small  particles,  the  curves  follow  a  X  law;  i.  e.  ,  the  Mie  theory 
reduces  to  the  Rayleigh  scattering  law  for  small  particles. 
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References  9  and  10  discuss  the  problem  of  scattering  by 
a  natural  haze.  Wavelengths  investigated  ranged  from  the  visible 
(0.  6l|i)  to  the  infrared  (11.  48}i).  Measurements  showed  that  the 
nuclei  which  caused  attenuation  of  visible  radiation  also  caused 
attenuation  at  infrared  wave  lengths.  According  to  Gebbie  most  of 
the  haze  droplets  have  radii  in  the  range  0.  ifi  to  0.  6|i;  he  assumes 

2 

the  distribution  function  of  the  radii  is  of  the  form  exp  -  k(l  -  r/r^) 
where  r^  is  the  most  probable  radius  (about  0.  4p).  Thus  the  size 
parameter,  a*  of  most  particles  in  a  natural  haze  for  visible  and 
longer  wavelengths  is  less  than  6.  Qualitatively  the  attenuation  by 
scattering  of  haz^^  can  be  explained  by  the  Mie  theory.  The  scat¬ 
tering  coefficient*  a  function  of  K,  has  a  peak  value  near  0.  61  p 
ajid  decreases  as  shown  in  Figure  22  with  increasing  wavelength; 
thus,  in  a  haze,  infrared  radiation  should  be  less  affected  than 
visible  radiation  and  long  wavelengths  should  be  scattered  according 
to  Rayhaigh's  law  -  which  corresponds  to  observations.  Exact 
calculations  are  impossible  because  of  the  lack  of  extensive  data 
on  the  size,  distribution,  and  characteristics  of  the  scattering 
particles. 

The  straight  line  in  Figure  23  represents  an  attempt  by 
Gibbons^ to  obtain  a  simple  approximate  expression  for  the 
scattering  coefficient  of  natural  haze.  (The  curve  must  lie  between 
the  curves  for  dry  salt  and  pure  water  as  those  curves  represent 
the  limiting  values  for  natural  droplets. )  Over  the  range  of  a 
which  is  appropriate  for  the  size  of  particles  involved,  the  approxi¬ 
mation  for  K  is  K  -  0.  9  Of  =  0,  6.  In  general,  the  scattering  coeffi¬ 
cient  CT  for  the  ith  type  particle,  of  which  there  are  n.  per  unit 
volume,  is 

a.  --  irr.^  N.K. 

i  lit 

K.  ~  a.  -  0.  67 

1  1 

a.  (2ifr./ K  -  0.  67/ X^) 

Considering  the  entire  haze: 
n  n 

Cr=2  ^i^'S  <2Trr./K)  -n(0.67/)^) 
i=l  i=l 


39 


From  this  relation  it  appears  that  a  relation  of  the  type 
(7=  c**'  where  -  1  <n <0 would\be  appropriate  for  v.'avelengths  less 
than  4p,  Gibbons  determined  tl^at  the  value  N  =  -0.  7  gave  a  good 
fit  to  Gebbie's  data  over  the  entire  range  of  measurements  -  even 
though  n  ehould  be  positive  for\arger  wavelengths.  Therefore 
scattering  by  natural  haze  can  be  assumed  to  follow  the  law: 

T  =  ‘  fcr  X.  from  visible  to  lOp,  and  follow 

the  Rayleigh  scattering  law  for  longer  wavelengths. 

Mie  Theory  can  also  give  a  qualitive  explanation  for  the 
attenuation  by  fog.  The  distribution  of  particle  sizes  in  a  fog  has 
a  form  similar  to  that  of  haze  particles,  except  that  there  are  many 
more  large  drops-radii  in  the  order  of  2  to  7fi.  This  means  that 
there  are  many  more  particles  with  size  parameter,  a,  relatively 
large.  Examining  the  curves  for  K,  one  sees  that  the  scattering, 
coefficient  will  be  appreciable  for  wavelengths  up  to  40  or  SOp. 
Droplets  in  clouds  and  in  rain  are  still  larger  so  that  Infrared 
energy  will  suffer  almost  neutral  extinction;  i.  e,  ,  for  infrared 
wavelengths  K  will  be  approximately  constant  (=  2),  so  that  scat¬ 
tering  will  be  complete  for  short  paths.  Figure  24  presents  the 
scattering  distribution  for  a  relatively  large  spherical  particle. 

The  predominance  of  scattering  in  the  forward  direction  is  called 
the  Mie  eff.ct.  It  has  beeu  found  that  this  effect  increases  with 
increasing  a* 

CLOUDS 

Clouds  play  an  important  role  in  the  determination  of  the 
background.  For  this  reason  they  are  considered  separately,  even 
though  they  are  only  special  cases  of  the  processes  already  con¬ 
sidered.  Low  and  middle  altitude  clouds  are  made  up  of  water 
droplets;  high  clouds  comast  of  ice  cryr^als.  Because  of  the  R 
relation  discovered  by  Junge'  '  the  most  important  nuclei  which 
act  as  condensation  centers  are  those  with  radii  0.  ip  to  1.  Ofx. 

These  nuclei  arise  from  two  main  sources.  Over  the  ocean  salt 
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particles  produced  by  the  wave  motion  of  the  sea,  and  particles 
produced  by  combustion  processes,  dust  storms,  and  volcanic 
action  on  land.  Those  nuclei  produced  on  land  are  initially  small, 
but  grow  by  coagulation  to  the  observed  ‘^ize.  In  the  troposphere, 
there  are  always  enough  nuclei  present  so  that  condensation  will 
occur  near  the  point  of  water  saturation'. •  The  growth  rate  of  a 
droplet  dR/dt  is  inversely  proportional  to  the  size  of  the  droplet. 
Thus,  even  though  there  may  be  a  cons id^'^r able  variance  in  the 
initial  size  of  the  nuclei  in  an  air  mass  which  rises  in  an  updraft, 
the  cumulus  cloud  which  forms  will  contain  a  narrower  spectrum 
of  droplet  sizes.  However,  as  the  cloud  continues  to  growi  there 
will  be  a  broadening  of  the  spectrum  due  to  coalescence.  Reference  8 
presents  several  theories  which  attempt  the  meci\anism  of  coales¬ 
cence,  but  the  process  is  not  well  defined.  Stratus  clouds,  formed 
when  a  layer  of  air  is  cooled  without  vertical  movement,  display 
the  same  type  of  spectrum  broadening  with  time.  Figure  25  pre¬ 
sents  droplet  spectra  for  various  types  of  low  and  middle  altitude 
clouds.  After  a  stratus  cloud  has  grown  to  a  few  hundred  meters 
in  thickness  it  will  start  to  drizzle.  Because  of  their  continual 
development,  it  in  impossible  to  consider  these  clouds  to  be  any 
given  state;  rather,  each  is  part  of  a  process.  For  this  reason  it 
is  impossible  to  obtain  an  exact  representation  in  the  model  of  these 
clovids.  They  will  be  considered  only  in  the  gross  sense  initially 
described;  their  mechanisms  of  scattering  and  absorbing  are  special 
cases  of  the  theories  discussed  above. 

In  contrast,  the  high  level  clouds  -  those  which  appear  above 
20,000  feet,  can  be  considered  to  be  in  a  relatively  stationary  state. 
These  clouds  are  classed  as  cirrus  -  which  includes  the  cirrostratus 
and  cirrocumulus,  nacreous,  and  noctilucent.  Only  a  very  limited 
amount  of  inforn  ation,  gathered  by  aircraft  penetrations,  is  avail¬ 
able  about  them,  fhe  accompanying  table  summarizes  the  char¬ 
acteristics  of  cirrus  clouds. 
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Cirrostratus 

Cirrocumulus 

Particle 

Ice  Crystal,  possibly  a  few 

Ice  crystal,  probably 

water  drops. 

some  water  drops. 

Shape 

Single  well  formed  column 

Bundles  of  incompletely 

or  twin  crystal. 

formed  columns. 

• 

Size 

100^  long,  40p  wide 

200-300^  long, 

50-100^  wide 

r' 

f 

Particle  Concen 

to®  -  10^ /m^ 

10^  -  10^ /m^ 

• 

tration 

Water  Content 

0.  01  g/m^ 

0.  1  -  0.  4  g/m^ 

r 

IR  Transmission 

Per  meter 

0.  992  -  0.  9992 

0.  53  -  0.98 

Per  100  meters 

0.46  -  0.92 

0  -  0.  09 

Per  1000  meters 

0.  0005  -  0.  46 

0 

r'  Clrrostratus  is  the  more  prevalent  type  of  cirrus  cloud, 

^  appearing  over  wide  areas  preceding  storms.  Cirrocumulus  may 

—  appear  with  cirrostratus  and  sometimes  with  thunderstorms.  The 

altitude  and  thicknes^i  of  these  clouds  is  a  function  of  season  and 
_  latitude.  For  the  purposes  of  the  model  average  values  will  be 

taken.  The  average  alt;.*ude  is  32,000  feet;  the  average  thickness 
is  2700  feet;  they  occur  about  15  percent  of  the  time. 

Nacreous  clouds  appear  in  the  arctic  region  at  altitudes  of 
21  to  30  km.  Their  frequency  of  occurrence  is  quite  low,  and  their 
infrared  transmission  probably  quite  high.  Noctilucent  clouds  appear 
only  rarely.  The  most  recent  data  indicates  that  they  appear  at  an 
average  neight  of  83  km  over  high  latitudes:  45^  -  71°  N.  Their 
transmission  is  probably  high.  In  the  model  these  two  types  of 
clouds  will  be  neglected;  cirrus  clouds  will  be  assumed  to  be  the 
only  high  clouds  of  significance. 
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OBSERVED  BACKGROUND  RADIATION 


The  following  table  summarizes  the  sources  which  contribute 
radiation  to  both  daytime  and  night  time  backgrounds: 


Radiation  Sources 


Daytime  Background 

1.  Scattered  radiation  (both  primary 
and  multiple),  from  the  sun. 

2.  Solar  radiation  reflected  directly 
to  an  observer  by  the  earth  or 
clouds. 

3.  Atmospheric  emission. and  earth 
emission  (for  simplicity  the 
earth  is  ai^sumed  uniform). 


Night  Time  Background 

1.  Radiation  from  celestial  bodies 
other  than  the  sun.  Assvimed 
negligible  for  simplicily. 

2.  Emission  by  clouds  and  the 
atmosphere. 

3.  Emission  from  the  earth 
which  is  assumed  to  be  a 
uniform  source;  i.  e.  , 
scattered  intense  sources  are 
assumed  to  average  out  over 
the  surface. 


Figure  20.  Atmospheric  Model 
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Figure  22.  Effective  Scattering  Ratio  to  Geometrical  Cross-Sectional  Area. 


Figure  23..  Scattering  Area  Coefficient  for  Several  Refraction  Indices. 
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Figure  24.  Scattering  Distribution  for  a  Large  Spherical  Particle 


Figure  25.  Dro;.4et  Spectra  for  Various  Types  of  Clouds. 
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EMPIRICAL  CONSIDERATIONS  FOR  A 
BRIGHTNESS  MODEL 


Delineated  in  the  previous  section  is  a  general  broad  outline 
of  the  severajl  processes  which  affect  the  infrared  background.  To 
formulate  a  completely  general  mathematical  model  even  for  the 
region  of  the  spectrim  from  1-10  microns  represents  a  highly  am¬ 
bitious  project.  To  be  accurate,  such  a  model  would  no  doubt  be 
complex.  If  it  is  to  be  more  than  a  concept,  and  thus  produce  real 
numbers,  severe  burdens  in  labor  and  time  may  be  imposed  in  its 
solution  for  specific  problems.  In  many  cases  this  is  not  justified 
by  the  requirements  and  a  reasonable  approximation  is  all  that  may 
be  required.  The  purpose  here  is  to  present  an  expression  of  pro¬ 
portionality  that  may  be  developed  into  a  reasonable  approximation 
of  the  sky  brightness  function.  Such  ah  expression,  because  of  its 
simplicity,  must  necessarily  ignore  the  many  causes  of  variations 
in  the  background  and  must  be  judged  on  the  merits  of  its  general 
or  average  accuracy.  While  the  data  taken  at  Key  West  has  been 
used  to  make  the  approximations,  the  formula  cannot  be  made  to 
match  at  every  point.  In  fact,  the  formula  can  only  be  made  to 
approximate  an  average  for  most  of  the  data.  It  must  be  remembered 
that  the  backgrounds  measured  at  Key  West  were  highly  complex  in 
their  natujre  and  are  not  themselves  representative  of  average  back- 
ground^j.  Measurements  fo  a  blue  cloudless  sky  are  required  to 
make  a  beginning  in  this  direction  and  it  is  hoped  that  the  propor¬ 
tionality  will  be  tested  against  such  ^rieasurements.  Investigation 
could  then  be  made  as  to  the  best  manner  of  adjustment  for  changing 
conditions. 

Upon  examination  of  background  measurements,  there  are 
some  trends  that  always  occur.  The  energy  received  by  a  measur¬ 
ing  radiometer  is  generated  by  two  sources,  the  earth  and  the  sun. 
The  energy  in  the  field  of  view  is  propagated  along  the  line  of  sight 
modified  by  the  complex  processes  of  primary  and  higher  order 

scattering,  absorption,  reflection,  refraction,  diffraction  and 

(17  18) 

emission.  In  general,  however,  spectral  measurements  '  ^  ' 


of  the  background  characteristically  show  two  sources  of  radiation, 

1)  nearly  ambient  temperature  radiation  (earth)  and  2)  high  tem¬ 
perature  radiation  (siin).  Some  of  the  low  temperature  radiation  is 
due  in  part  to  atmospheric  constituents  absorbing  high  temperature 
solar  radiation  and  in  turn  re -emitting  as  a  low  temperature  source. 

Measurements  made  in  terms  of  elevation  and  azimuth  of  the 

observed  point  also  show  certain  general  characteristics.  It  is 

noticed  that  the  radiance  emitted  by  clear  sky  is  higher  near  the 

horizon  than  anywhere  else  and  decreases  with  elevation.  Again, 

when  one  measures  the  sky  towards  the  sun,  the  radiation  levels 

(19) 

increase  rapidly  closer  to  the  sun.  Bennett,  Bennett  and  Nagle'  ' 
show  that  the  measured  radiance  at  the  horizon  is  closely  that  from 
a  blackbody  at  approximately  ambient  temperature. 

Depending  upon  the  position  of  the  sun^  there  would  be  addi¬ 
tionally  an  increase  at  shorter  wavelengths  (i-Sp)  due  to  increased 
back -scattering  from  the  denser  air  mass  as  one  measures  closer 
to  the  horizon.  This  occurs  when  the  sun  is  behind  the  observer 
so  that  the  observed  point  is  at  a  large  angle  with  the  sun  which  is 
primarily  in  early  morning  and  late  afternoon  hours.  A  further 
condition  exists  due  to  polarization.  A  minimum  is  expected 
at  90^  from  the  sun  but  this  is  usually  found  to  be  shifted  towards 
the  zenith.  Similar  trends  can  be  noted  in  the  Polar  Plots  presented 
in  the  first  part  of  this  report,  even  though  the  sky  is  overcast  or. 
at  best,  highly  clouded,  and  the  trends  can.  be  due  to  other  causer. 

A  plot  of  elevation  with  constant  azimuth  for  the  overcast  sky  shows 
the  radiation  levels  near  the  horizon  to  be  c..ose  to  the  blackbody 
levels  at  ambient  temperature,  but  with  increasing  levels  encountered 
with  elevation  angle  in  some  cases  due  no  doubt  to  change  in  tempera¬ 
ture  of  the  cloud  cover,  or  differences  in  ambient  temperature  at 
higher  altitudes.  This  would  be  apparent  in  total  radiation  measure¬ 
ments  that  cover  a  broad  spectral  range  as  thcfe  do,  0-l6p.  For 
the  4.  25  micron  region,  if  the  data  shown  in  Figure  4  is  smoothed 
to  remove  some  of  the  detail,  it  would  appear  as  in  Figure  26.  No 
minimum  is  observed  due  to  polarization  but  the  trend  of  decrease 
with  elevation  from  the  horizon  and  increase  toward  the  sun  can  be 
noted. 
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In  Figure  5,  which  is  for  blue  sky  of  June  23rd,  a  minimxim 
is  noted  approximately  60^  from  the  sun  and  the  usual  decrease 
from  horizop  and  increase  toward  sun  as  well.  It  will  be  repeated 
here  that  these  affects  probably  would  not  be  so  apparent  except 
for  the  fact  ^at  the  filter  used  to  make  the  measurements  has  a 
bandwidth  from  3.75  to  4.75  microns  and  consequently  admits 
radiation  outside  of  the  CO2  absorption  rej;ion.  These  effects  notdd 
are  probably  due  more  to  this  radiation  than  to  CO2  emission. 

To  attempt  a  computation  which  would  construct  these  polar 
diagrams  synthetically,  two  general  shaping  functions  have  been 
chosen,  one  to  represent  the  family  of  curves  originating  from  tie 
sun  position  and  the  other  to  represent  the  family  of  curves  origi¬ 
nating  at  the  horizon.  The.  sum  of  these  two  functions  would  then 
represent  a  proportionality  which,  when  used  with  suitable  constants, 
might  approach  the  actual  radiation  measurements  of  the  sky.  Sim¬ 
plicity  must  also  be  the  keynote  for  such  functions  or  their  useful¬ 
ness  is  completely  nulled.  Since  they  are  approximations  only,  a 
great  deal  of  time  ior  computing  should  not  be  required. 

For  severc.1  pieces  of  data  at  Key  West,  the  family  of  circles 
originating  with  the  horizon  can  be  approximated  by 


Where  G  is  the  elevation  angle,  at  any  azimuth,  measured 
from  the  horizon. 

For  the  same  pieces  of  data,  the  family  of  circles  originating 
with  the  sun  can  be  approximated  by 


1 


(2) 


Where  0  is  the  angle  between  the  sun  and  the  observed  point 
measured  from  the  observer’s  station.  Figure  27  is  an  illustration 
in  point  for  the  0«l6  micron  region  for  two  days.  These  functions 
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were  used  to  compute  the  same  points  mt^asured  and  published  in 
table  form  in  Repc.L^t  No.  3  for  this  contract.  '  '  The  original  data, 
the  computation  results,  and  the  percentage  differencec  are  pre¬ 
sented  in  Tables  III,  IV,  and  V,  and  polar  plots  have  been  made 
from  the  computations  to  compare  with  the  data  plots.  Figures  Z8, 
29,  and  30  are  the  computed  diagrams  to  be  compared  with  Figures 
2,  3,  and  5,  respectively. 

The  equation  was  written  as  follows: 


(3) 


where  N,  is  the  radiance  of  the  background,  N  's  the  radiance 
due  to  ambient  temperature  and  is  the  radiance  due  to  the  sun. 
Since  it  is  not  possible  at  this  time  to  present  a  finished 
formula  from  the  Key  West  data,  but  only  to  establish  a  proportion¬ 
ality,  numbers  were  chosen  for  N  and  N  to  fit  the  data.  N 

as  a 

was  taken  as  the  blackbody  radiance  for  ambient  temperature  for 
the  wavelength  region  under  consideration  and  this  works  out  very 
well.  In  a  sky  completely  overcast  with  clouds,  it  would  be  almost 
impossible,  without  either  measurement  or  extended  analysis,  to 
choose  a  value  for  N  .  The  0-16  micron  region  contains  almost 
all  of  the  radiation  for  blackbody  of  5500°C  which  is  the  equivalent 
for  the  sun.  A  great  deal  of  this  radiation  would  be  lost  due  to 
backscattering  ^nd  reflection  from  the  top  of  the  cloud  cover.  In 
turn,  a  large  amount  would  be  absorbed  and  re  - radiated  with  ambient 
temperature  distribution.  There  would  be  still  further  attenuation 
due  to  high  relative  humidity  at  ground  level.  It  was  evident  that 

visible  illumination  v/as  extremely  low  on  June  16  and  17.  was 

4  2  ^ 

taken  as  4.45  x  10  fjiw/cm  .  5^.eradian  for  the  0-16  micron  region 

and  178  jjLw/cm^.  steradian-  micron  for  the  4.  25  micron  region. 
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To  N  was  added  12  percent  of  N  as  a  coefficient  for  f(0) 
because  higher  order  scattering  causes  bright  sky  along  the  hori¬ 
zon  resulting  in  uniform  intensity  with  respect  to  azimuth, 
and  for  the  present  is  assumed  to  be  a  function  of  elevation  from 
the  horizon.  These  values  were  chosen  to  suit  the  0-l6  micron 
data  for  June  l6th.  Table  III  gives  the  observed  values  made  dur¬ 
ing  the  field  trip,  and  below  these  values,  computed  with  equation  (3), 
The  percentage  error  has  been  obtained  simply  by  taking  the  dif¬ 
ference  between  observed  and  computed  values  and  dividing  by  the 
observed  value.  Whether  the  error  is  plus  or  minus  indicates  the 
computation  was  high  or  low.  For  the  25  observed  points,  there 
are  14  minus  values  and  11  plus  values  with  a  median  around  -2 
percent.  The  highest  deviation  was  plus  22  percent  and  the  lowest 
was  minus  35  percent.  Figure  28  is  a  polar  plot  of  the  computation 
and  is  to  be  compared  with  the  plot  of  the  actual  data.  Figure  2. 

The  same  values  for  N  and  N  were  used  to  repeat  this 

cL  S 

procedure  for  the  0-16 micron  region  for  the  following  day.  The 
results  are  shown  in  Table  IV  and  the  computation  is  plotted  in 
Figure  29  to  be  compared  with  the  actual  data  plot  of  Figure  3. 

Here  the  computation  is  a  little  on  the  high  side  with  a  median  of 
plus  14  percent.  A  very  small  downward  revision  of  N  would 
have  corrected  this  since  June  17th  was  darker  than  June  l6th. 

Table  V  and  Figure  30  are  the  results  for  the  4.  25  region 
giving  a  median  error  of  plus  1-1/2  percent  with  more  than  1/2 
the  values  being  on  the  plus  side.  These  are  to  be  compared  with 
the  actual  data  plot  of  Figure  5.  A  computation  was  made  but  not 
plotted  for  the  4.  25  region  for  June  17th  with  the  results  in  Table 
VL  These  can  be  compared  to  the  actual  data  plot  of  Figure  4. 

As  in  Table  IV  for  0-16  microns,  the  values  are  all  on  the  plus  side 
with  a  mean' of  approximately  plus  38  percent. 

It^ppears  that  a  proportionality  exists  which  may  be  simple 
to  compute  and  accurate  enough  for  many  engineering  problems. 

In  a  cloudless  and  clear  sky,  the  radiance  levels  near  the  sun  have 
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a  much  steeper  gradient  than  in  a  hazy  sky  where  the  increased 
scattering  would  tend  to  broaden  the  gradient.  For  the  same  Treasons, 
the  apparent  emissivity  of  the  atmosphere  would  tend  to  flatten  as 
the  air  became  filled  with  more  scattering  particles. 

If  f(3t  ip  written 


(4) 


and  plotted,  as  in  Figure  31,  for  different  values  of  jS,  the  form  is 

(191 

similar  to  that  given  by  the  expression  for  the  emissivity  of  the  sky' 
when  ^  >  1.  As  ^  increases,  the  resultant  curve  approximates 
emission  curves  for  the  sky  as  the  sky  becomes  more  blue  with 
less  scattering  in  the  near  infrared  as  a  result. 

In  the  same  way,  f(^)  can  be  written 

(5) 


with  the  result  that,  as  a.  increases,  the  steepness  of  the  gradient 
increases  about  0  =  0^,  as  shown  in  Figure  32.  and  a»  therefore 
can  be  taken  as  "indicators"  of  the  condition  of  the  sky,  becoming 
larger  as  the  sky  becomes  clearer  with  less  scatter.  When  there 
are  clouds  present  in  the  sky,  polar  diagrams  would,  in  certain 
wavelength  regions  more  than  others,  produce  higher  values  than 
the  blue  sky,  and  produce  data  that  was  very  erratic  in  form  be¬ 
cause  it  would  represent  a  large  amount  of  detail.  Under  the  over¬ 
cast  conditions  at  Key  West,  plots  of  energy  levels  with  respect  to 
0  but  at  constant  azimuth  show  a  convex  appearance  that  is  no 
doubt  due  to  heating  of  the  cloud  cover  or  a  gradient  in  air  tempera¬ 
ture  with  altitude.  If  this  is  a  general  condition  for  certain  types 
of  sky,  then  fractional  values  of  j?  will  approximate  the  condition 
a.i  seen  in  Figure  31.  The  detail  of  the  cloud  structure  would  disap¬ 
pear,  but  an  average  value  would  result. 
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No  doubt  if  a  and  were  to  be  evaluated,  it  would  be  seen 
that  they  are  really  complex  functions.  Certainly  they  are  a  func¬ 
tion  of  wavelength,  and  measurements  on  blue  uky  would  determine 
values  for  specific  wavelength  regions.  All  cf  the  computations 
given  here  were  for  ^-1.5  and  a  =  1.  4,  or  approximately  equal. 

It  may  be  that  one  index  of  the  sky  is  all  that  is  necessary  to  give 
a  good  approximation.  Measurements  of  a  and  /?  would  also  pro¬ 
duce  the  proper  values  for  N  which  would  be  more  compatible 
with  expected  values,  especially  in  a  clear  sky. 

The  computations  weye  not  very  difficult.  The  only  one  which 
appears  troublesome  is  the  determination  of  This  is  a  function 
of  the  elevation  and  azimuth  of  the  sun,  and  the  elevation  and  azi¬ 
muth  of  the  observed  point,  and  is  expressed  as  follows: 

cos  ^  =  sin  w  sin  9+  cos  w  cos  0  cos  p  (6) 


where 


p  =  the  difference  in  azimuth  of  sun  and  observed  point 
w  =  sun  elevation 
9  =  observed  point  elevation 

Most  values  of  0  used  in  the  computations  were  read  by  using 
dividers  on  a  ball  having  latitude  and  longitude  lines  every  10  de¬ 
grees.  The  dividers  measure  the  chord  distance  between  two 
points  and  this  is  used  on  any  longitude  (great  circle)  line  to  read 
degrees,  0  has  been  read  to  one  degree  accuracy  in  this  fashion 
on  a  9 -inch  diameter  ball  which  is  accurate  enough  for  ^  greater 
than  30^.  When  0  is  smaller,  or  changes  slowly,  as  in  the  case 
for  target  azimuth  planes  at  small  angles  to  the  sun  azimuth  plane, 
<f>  must  be  computed  with  equation  (6).  The  values  for  the  sines  of 
9  and  0  in  equation  (3)  were  to  3  decimals  and  the  calculations  were 
performed  on  a  slide  rule. 
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minus  27  plus  median  approximately  +14% 
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32  plus  values  -  3  minus  values  -  median  approximately  +38% 


Overcast  Sky,  June, .17,  1959,  4.  25  Microns 
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Figure  30.  Computed  4.  25  Micron  Region  for  June  23rd. 


Fig'iire  31.  Effect  of  on  f(G) 
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Figure  32.  Effect  of  a  on  f(^) 
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ESTIMATED  EARTH  AMD  SKY  BACKGROUND 
NOISE  INTERFERENCE  IN 


INFRARED  DEVICES 


INTRODUCTION 

Backgrounds  are  a  con:iplex  structure  which  reflect,  scatter, 
and  absorb  incident  electromagnetic  radiant  energy  and,  further, 
reemit  some  of  this  radiant  energy  due  to  their  own  thermal  agita¬ 
tion  and  chemical  composition.  This  gaseous  composition  of  the 
atmosphere  varies  further  with  time,  altitude  and  location  on  the 
earth.  However,  if  one  disregards  the  source  from  which  the  elec¬ 
tromagnetic  energy  may  have  arisen  and  the  corresponding  operations 
through  which  it  may  have  been  subjected,  some  meaningful  char¬ 
acteristics  of  backgrounds  may  be  noted  based  upon  strictly  empirical 
observations.  The  background  is  there  and  possesses  certain  dis¬ 
tinguishable  features  character  ized  by  the  spatial  structure  and 
spectral  distribution  of  the  radiation.  Some  of  these  salient  char¬ 
acteristics  are  briefly  noted  below  in  Radiometric  Observations. 

In  order  to  more  fully  and  conveniently  apply  empirical  gra¬ 
dient  data  regarding  the  background,  it  has  proven  useful  to  develop 
a  hypothetical  one  dimensional  model  of  background  interference 
described  belov/  in  One  Dimensional  Background  Noise  Spectra. 

This  model  has  proven  especially  valuable  in  indicating  expected 
interference  distribution  in  regions  of  the  spectral  or  spatial  dis¬ 
tribution  where  measurements  have  not  as  yet  been  made  in  great 
detail--for  example,  at  high  space  frequencies.  This  heuristic 
model  is  developed  below  and  indicates,  for  example,  the  impor¬ 
tance  of  contrast  in  determining  background  interference,  and 
further,  indicates  the  expected  space  frequency  accelerated  drop 
off  at  high  frequencies,  and  the  expected  levelling  off  at  low  space 
frequencies  (it  cannot  continue  to  be  1/f^  to  zero  space  frequency). 
This  spatial  structure  of  the  power  density  function  will  also  be 
shown  to  vary  with  altitude  or  distance  to  interfering  backgrounds 
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(clouds)  so  that  satellites  viewing  the  earth  will  prei^u-rnably  see  a 
different  spatial  structure  than  we  are  accustomed  to,  i.  e.  scaled 
down  but  with  identical  contrasts.  This  expected  variation  in 
structure  with  angle  and  altitude  is  described  below  and  will  be 
corroborated  in  the  near  future  by  a  detailed  examination  and  corre¬ 
lation  of  TIR^  satellite  pictures  of  the  earth  within  the  limits  of 
the  resolving  power  of  that  instrument. 

The  one  dimensional  Wiener  spectrum  of  background  noise 
has  only  limited  usefulness  in  the  calculation  of  noise  interference 
levels  when  used  with  two  dimensional  field  stops  and  reticles. 
Unfortunately,  most  background  gradient  structure  is  obtained  in 
this  manner  today  despite  methods  proposed  by  R.  Clark  Jones 
several  years  ago  to  obtain  the  two  dimensional  spectra  directly.  '  ' 

A  method  of  reducing  the  one  dimensional  data  to  a  two  dimensional 
Wiener  spect’.a  is  described  in  Two  Dim  Vinal  Backgr Noise 
Spectra. 

Finally,  one  may  apply  directly  the  two  c  .  /  noise 

spectra  against  the  field  stop  and  reticle  to  obtain  i  estimate  of 

(25) 

expected  background  interference. '  '  This  method  is  illustrated 

briefly  in  the  section  on  Reticle -Aperture  Space  Filtering. 

RADIOMETRIC  OBSERVATIONS 

Figure  33  illustrates  some  typical  ambient  radiometric  meas¬ 
ures  of  background  intensity  as  a  function  of  wavelength.  Back¬ 
ground  interference  existing  primarily  in  wavelengths  shorter  than 
3p  occur  as  a  result  of  diffusely  reflected  and  scattered  sunlight. 
Atmospheric  absorption  bands  due  to  CO^  and  H^O  are  evident. 

Re -emission  of  solar  energy  in  the  absorption  bands  occurs  at  much 
lower  temperatures  than  the  apparent  5900°K  sun  and  does  not  fill 
the  absorption  band  corresponding  to  the  decrease  in  solar  intensity. 

In  the  thermal  emission  band,  beyond  4(i,  the  average  spectral 
distribution  follows  that  of  an  apparent  300°K  blackbody.  In  this 
spectral  band,  atmospheric  absorption  and  re -radiation  occur  at 


closely  the  same  apparent  temperature  and  a  corresponding  smooth 
emission  curve  results.  Beyond  8p.,  however,  little  atmospheric 
emission  occurs  because  the  longer  wavelength  involved  is  ifi- 
cantly  larger  than  the  air  molecules  themselves  (little  atmospheric 
scattering  is  experienced  at  these  longer  wavelengths).  Thus,  only 
emission  from  clouds  beyond  8fi  permits  a  300^K  blackbody  approxi¬ 
mation  of  background  interference  when  looking  up.  When  looking 
down  at  the  earth  between  8-13p.  there  would  presumably  be  less 
contrast. 

It  will  be  developed  shortly,  that  absolute  background  contrasts 
are  the  greatest  single  factor  in  developing  interference  levels  in 
various  spectral  regions.  We  are  familiar  with  the  existence  of 
relative  contrasts  in  the  visible  spectrum  and  the  enhancement  of 
these  contrasts  (at  least  for  clouds)  by  the  use  of  near  infrared 
film  neiJV  Ip.  Actual  contrast  measurements  betv'een  clouds  and 
sky  are  shown  to  some  extent  in  Figure  33.  These  contrasts  are 
moderate  towards  the  blue  end  of  the  visible  spectrum,  and  then 
increase  to  ratios  of  10,  or  100,  or  a  1000  to  1  in  the  I  -  3p  band. 

It  is  these  large  contrasts  that  make  daylight  background  interference 
in  the  1  -  3p  band  such  a  formidable  problem.  Contrasts  in  the 
longer  wavelength  regions  between  3.  5p  and  8p  are  small  -  an  ideal 
region  for  obtaining  minimum  background  interference.  From  8  -  13|i 
large  contrasts  exist  again  between  cloud  and  clear  sky  indicating 
another  regicn  of  large  background  interference. 

ONE  DIMENSIONAL  BACKGROUND 
NOISE  SPECTRA 

The  variations  in  intensity  or  detail  of  light  distribution  which 
arise  L'om  the  background  in  space,  has  been  described  by  R.  Clark 
Jones  as  space  power  -  or  Wiener  -  spectra  with  ordinates  of  normal¬ 
ized  space  frequency.  This  method  of  describing  actual  backgrounds 
has  been  used  with  success  in  gradient  meter  measurements  similar 
to  the  one  conducted  by  the  The  Ramo-Wooldridge  Corporation  under 
Contract  AF  19(604)  3473  from  AFCRC  and  earlier  under  contract 
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AF  33(600)  30489.^^*'^^  These  spectra  describe  the  background  as 
varying  approx: niiately  inversely  to  some  power  of  the  space  fre¬ 
quency  -  at  l^ast  at  tl-e  higher  space  frequencies  that  hav^  been 
measured.  But  this  method  of  measurement  gives  only  a  partial 
description  of  the  background  since  dll  phase  information  is  lost; 
i,  e.  ,  an  inverpe  ypace  frequency  spectra  itself  may  arise  Irom 
different  background  configurai:icn8.  Fv  rthermore,  as  a  practical 
matter,  if  it  is  desirable  to  distinguish  cr  highli^jht  a  particular 
size  or  shape  object  against  the  background,  it  may  require  differ¬ 
ent  optimum  filtering  techniques  -  such  as  a  discrete  or  non-linear 
space  filter. 

The  use  of  unrestricted  theoretical  models,  nowc/er.  may 
suggest  alternate  study,  measurement  or  testing  methodj.  For 
example,  the  model  presented  describes  the  nature  of  ihe  rransitioiA 
from  the  ambient  intensity  level  to  the  lowest  space  frequency 
spectra  whicn  can  practically  be  measured.  The  model  also  per¬ 
mits  an  estimate  of  higher  probability  distribution  of  the  background 
and  permits  alternate  improved  filtering  techniques.  Finally,  the 
use  of  a  theoretical  model  may  permit  a  systematic  method  of 
normalizing  or  ciandar^lizing  background  space  filtering  techniques 
and  their  methods  of  tetit. 

It  will  be  assumed  in  the  discussion  which  follows  that,  by 
and  large,  the  greatest  single  source  of  background  noise  occurs 
from  cloud  structure  -  either  against  the  sky  v/hen  looking  up  - 
or  against  the  earth  when  looking  down.  A  Poisson  distribution  of 
cloud  sizes  will  be  assumed  with  a  50  percent  probability  of  cloud 
occur  ranee  at  any  point  around  the  world.  This  assumption  com¬ 
pares  with  an  actual  average  52  percent  cloud  cover  around  the 

/  5  \ 

world,  according  to  C.  E.  P.  Brooks.  '  ' 

A  series  of  pulses  of  random  spacing  and  random  width  are 
used,  all  with  a  constant  amplitude  or  height.  These  pulses  then 
represent  clouds  against  a  biuu  sky  (or  earth)  when  scanned  with 
an  infinitesimally  sir  all  field  of  view  in  one  dimension.  Only  two 
amplitude  levels  are  assumed  to  exist.  The  high  level  represents 
reflected  sunlight  (or  ambient  temperature  of  the  cloud)  and  the 


lower  level  reprcEente  scattered  sunlight  or  earth  {or  apparent 
ambient  temperature  of  the  blue  sky).  The  distribution  of  pulse 
lengths  or  cloud  widths  P(d/  is  assumed  to  be  a  Poisson  distribution 
peaking  at  lengths  4  nau  tical  miles,  shown  in  Figure  34  and 
decreasing  in  probability  in  either  direction  of  length  or  width 
from  this  value.  The  use  of  a  particular  probability  of  pulse  dura¬ 
tion,  allows  us  to  derive  spectra  for  satellites  since  the  cloud 
angular  subtense  will  be  proportionally  reduced  with  altitude.  A 
more  exact  distribution  of  cloud  lengths  will  be  determined  in  the 
near  future;  hopefully  at  Ramo -Wooldridge ,  with  a  detailed  exami¬ 
nation  and  compilation  of  TIROS  pictures. 

Most  modern  probability  and  statistic  textbooks  describe  the 
mathematical  relations  and  functions  necessary  to  define  the  spec¬ 
trum  of  a  random  or  stochastic  process.  The  results  of  these 
derivations  only  are  presented  in  this  section  for  the  example  chosen 
to  represent  the  background. 

The  model  selected  to  describe  the  background  is  illustrated 
below.  The  length  of  pulse  or  spacing  is  distributed  in  space  x 
according  to  a  probability  function  P(a),  with  a  height  of  magnitude 
either  +1  or  -1.  In  practice,  the  difference  between  the  two  levels 
represents  contrast.  The  resulting  spectra  may  be  multiplied  by 
the  contrast  intensity  in  any  spectral  region  as  a  first  approximation 
to  the  gradient  background  interference. 


+  1 
o 

-  1 

Figure  35 


The  model  at  this  tixx.e  neglects  the  ambient  or  dc  level  of 
the  intensity  which  is  always  some  value  greater  than  zero.  How¬ 
ever,  the  average  level  may  be  readily  added  -  ard  this  will  be 
done  later  -  by  employing  the  superposition  concept  to  take  care 
of  the  ambient  level.  The  resulting  ambient  term  will  always 
appear  as  a  6  function  at  the  origin  with  an  integral  equal  to  the 
ambient  level. 

If  the  functions  0(f)  and  are  defined  as  follows: 


0(f)^  + 


P(a)e 


-j  Zirf 


then  the  normalized  power  spectral  distribution  is  given  by  H.  M. 
Foley  in  "An  Investigation  in  the  General  Theory  of  Pressure 
Broadening  and  an  Experimental  Study  of  Pressure  Effects  in  the 
14p  Band  of  Hydrogen  Cyanide.  "  PhD  Thesis,  University  of 
Michigan,  1942, 


S(f)  = 


1 


where  the  average  length  of  pulse  a  is 


nOO 

aP(a)d{a)  (3) 

The  space  frequency  ordinate  against  which  the  power  spectra 
will  be  plotted  is  f^. 

If  the  cloud  length  is  some  arbitrary  amount,  then  the  angular 
subtense  as  observed  from  afar  is  inversely  proportional  to  the 
distance  from  that  cloud.  For  convenience,  it  is  assumed  that  ail 
clouds  are  grouped  at  30,000  feet.  The  average  angular  subtense 
is,  therefore,  4/r  where  r  is  the  slant  range  in  nautical  miles  from 
oDserving  station  to  the  clouds.  The  probability  of  a  particular 
pulse  length  a  is  then 

P(a)  =  (4) 
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Solving  for  0(1)  and  ^f)  by  iBubstitutmg  Eqaativ^n  (4)  in 
Equation  (1),  we  have 


0(tl  4  (f)  * 


(r/4  +  j 


and  the  resub'-ing power  epectrvun  from  Equation  (2)  is 


S(y  =  ri 


o.^.r 


/ 


(£■; 


(6) 


since  “a-*  8/r. 

This  power  spectrum  ie  plotted  in  Figure  36  and  shows  a 
power  spectrum  "break  point"  which  varies  with  range  or  distance 
to  the  clouds.  Tne  dc  or  ambient  radiometer  results  in  a  dqlta 
iuncticr;  at  the  origin.  The  ambient  intensity  remains  invariant 
with  altitude,  or  aspect,  or  space  frequency  bandwidth  of  field  of 
view. 


It  ip  helpful  compare  the  ai/ove  model  with  some  actual 
sky  or  backgrcuryj  observations  and  measurements.  Several  aspects 
will  be  cone^idered;  (1),  the  relation  of  the  average  or  ambient  level 
with  respect  tu  the  contract  in  intensity  for  clouds  and  sky  back¬ 
ground  in  thf  near  anH  middle,  infrared;  (2),  the  space  frequency 
from  which  one  can  assume  that  the  .t/f  '  power  spectrum  intensity 
dependence  holds;  and  (3). the  effer.i  of  the  Rhsrpness  of  the  step 
changes  of  intensity  on  the  poweo  spectrum. 

In  the  n***ai  i;'irared  the  variatio  ,  in  intensity  between  a  white 
cloud  and  adjacent  blue  ?V.v  is  very  pronounced.  There  is  a  change 
of  an  order  of  magnitud.i  of  10  to  iOf'O  times  the  intensity  level. 

This  sharper ’’ng  of  contrast  in  the  near  infrared  has  been  frequently 
measured  and  is  borne  out  by  simple  observations  of  heightened 
cloud  contrast  in  photographs  taken  in  the  near  infrared.  The  varia¬ 
tion  of  intensity  in  the  mathematical  model  must  therefore  be  as¬ 
sumed  tc  be  of  the  same  order  of  magnitude  as  the  existing  ambient 
level  determined  from  the  solar  constaur  and  albedo  of  the  Earth. 


In  the  mi'v’Ue  infrared  from  3.  5  -  8|ji  where  the  background 
energy  primarily  arises  from  the  thermal  emission  of  cl^juds  or 
background,  however,  the  contrast  is  ,.iot  as  pronounced.  For 
example,  the' apparent  change  in  radiation  from  the  blue  sky  to  a 
cloud  is  of  the  order  of  20  percent  greater  than  the  blue  sky  thermal 
radiation.  The  ratio  of  the  intensity  level  change  for  the  variation 
from  the  mean  level  between  cloud  to  sky  is  therefore  of  the  order 
of  1/10  times  the  average  or  ambient  level.  In  the  8  -  13p  band, 
the  contra  re  are  larger  than  the  4  -  5p  band,  but  not  as  great  as 
the  1.  3  -  2.  7p  band.  Contrast  levels  equal  to  50  percent  of  the 
ambient  may  be  assumed  in  this  band.  The  predicted  background 
intensity  in  the  thermal  emission  band  beyond  3p  therefore  is  simply 
taken  as  the  product  of  the  proper  curve  of  Figure  36  times  one- 
half  the  contrast  level  or  0.  1  and  the  total  300°K  blackbody  radiation 
within  ^  particular  spectral  band.  The  procedure  for  estimating 
background  interference  in  the  shorter  wavelengths  during  daylight 
illumination  is  similar,  but  the  contrast  levels  are  assumed  to  be 
100  peresnt,  or  0.  50  for  1/2  the  contrast  level,  and  5  x  10  of  a 
5900^K  blackbody  is  used  for  solar  illumino.tion. 

Figures  37and  38  illustrate  examples  of  the  predicted  one 
dimensional  background  interference  as  well  as  actual  gradient 
ineasurements  obtained  in  the  field.  Reasonable  agreement  in  inter¬ 
ference  levels  is  evident. 

The  step  change  in  intensity  ac-sumed  in  the  calculations 
moreover,  is  net  realized  in  practice  -  especially  in  thermal 
emission  band.  This  fuzziness  of  the  edges  of  clouds  (or  edge 
softening)  will  further  attenuate  the  higher  space  frequencies.  For 
example,  if  a  30  foot  distance  at  the  edge  of  the  cloud  is  lequired 
to  go  from  one  extreme  intensity  peak  to  63  percent  of  the  value  of 
the  opposite  peak,  then  an  additional  break  point  in  the  power 
spectru^  N(f  )  will  occur  beyond  which  the  attenuation  is  1/f^. 

X  X 

This  cloud  edge  transient  .nay  commonly  occur  in  the  nea»  infrared. 

In  the  thermal  band,  however,  the  transient  length  is  probably  longer 
to  obtain  63  percent  of  the  peak  change  from  one  level  to  another  which 


will  thereby  shift  the  second  break  point  to  an  even  lower  space 
frequency.  There  has  been,  however,  insufficient  data  accumulated 
to  date  to  conclusively  verify  this  increased  attenuation,  but  it  is 
strongly  indicated  in  the  me  a,  ired  data  at  high  rpace  frequencies. 

A  clear  indication  of  any  such  trend  will  be  obtained  when  back¬ 
ground  data  is  obtained  for  space  frequencies  substantially  greater 
than  those  being  measured  today. 


TWO  DIMENSIONSAX.  BACKGROUND 
NOISE  SPECTRA 


The  background  noise  hac  been  considered  oniy  for  a  one 
dimensional  power  spectrum,  '^^2  conveision  from  one 

dimension  to  a  two  dimensional  noise  spectrum  is  considered  here 
as  well  as  methods  for  measuring  the  two  dimensional  spectra 
directly,  nN/^  +  fn 

An  aperture  with  sides  Ax  and  Ay  has  a  two  dimensional  filter 

function  r(f  f  )  defined  by 
X  y 


sin  irf  Ax 

_ X 

AX 


sin  Trfy.  Ay 
■nrfy  Ay 


If  Ax  and  Ay  are  approximately  equal  to  each  other,  one 
obtains  the  one  dimensional  noise  spectra  considered  earlier  after 
compensation  for  the  aperture  in  the  scan  direction.  This  correc¬ 
tion  factor  Is  the  reciprocal  of  the  apertv^re  function  in  the  scan 
direction  (if  the  filter  functionan  the  normal  direction  to  scan  is 
invariant)  i.  e.  ,  for  scans  in  x  direction  with  a  squars^  field  at 
space  freauency  fx. 

Aperture  Correction  =  - — 

sin  ir  f^  Ax 

This  aperture  when  convolved  with  the  background  noise, 
hov/ever,  is  then  represented  as  the  noise  existing  at  a  particular 
space  frequency  f  •  Actually,  the  background  noise  is  for  a  con- 
stant  f^  but  includes  a  wide  space  frequency  bandwidtn  in  the  normal 
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or  £  direction,  most  of  which  exists  at  distances  substantially 
y 

greater  than  from  the  origin  where  relatively  little  background 

occurs.  Atwodinxensional  spectra  derived  from  a  one  dimensional 

scan  ;zseumes  an  isotropic  behavior,  i.  e.  ,  that  the  background 

noise  level  is  only  a  function  of  \jf^  +  f  To  correct  the  one 

X  y 

dimensional  noise  spectra  for  two  dimensions*  one  must  then  divide 
the  resulting  output  by  the  effective  bandwidth  in  the  normal  direc¬ 
tion  to  the  scan. 

If  the  effective  bandwidth  is  considered  Vv)  lie  betw'een  the  1/2 
power  points,  then  in  the  presence  of  background  noise  varying 
inversely  with  space  frequency,  the  effective  bandw'idth  A  y  at  a 
one  dimensional  space  frequency  f^  is  A  y  =  l/2f^. 

At  the  low  space  frequency  end,  however,  the  background 
noise  spectra  is  constant  or  ’’white.  ”  In  this  region,  the  effective 
bandwidth  A  y  will  be  constant  and  approximately  equal  to  l/2^/T*  fb 
where  the  break  point  fb  is  the  space  frequency  sit  which  the  trainsi- 
tion  between  white  and  inverse  space  frequency  dependence  occurs. 
The  resulting  two  dimensional  noise  spectra. is  illustrated  in 
Figure  39. 

The  two  dimensional  characteristic  noise  spectra  of  back¬ 
ground  may  be  derived  directly  with  a  slit  scan  if  the  effective 
space  frequency  bandwidth  normal  to  the  scan  direction  is  small  in 

corx'iparjSon  with  the  space  frequency  in  the  scan  direction.  This 

(24) 

was  discussed  by  R.  Cl?rk  Jones  in  195  5'  '  for  regions  of  space 

frequency  f^  which  are  significantly  higher  than  1/2  Ay.  In  such 
circumstances,  the  effective  bandwidth  A  .y  in  the  direction  normal 
to  the  scan  direction  is  a  constant  and  the  resulting  two  dimensional 
noise  spectra  may  be  represented  as  occuring  at  space  frequencies 
equal  tc  _  ,  i.  e.  ,  I  ’  (~  if  L 
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RETICLE  -  APERTURTi:  SPACE 
FILTERING 


There  are  two  convenient  and  alternative  methods  of  analysis 
which  permit  a  study  of  the  space  discriminating  ability  of  optical 
systems --space  and  space  frequency.  The  detail  of  light  distribu¬ 
tion  across  the  object  or  image,  whether  from  the  background  or 
target,  can  be  expresses,  mathematically  as  a  function  of  its  inten¬ 
sity  and  the  linear  position  in  that  plane,  such  as  x  and  y.  It  may 
?lso  be  expressed  mathematically  as  a  function  of  inverse  space 
or  1/x  and  1/y,  in  which  case  it  is  analogous  to  frequency  in  the 
true  sense,  and  is  called  space  frequency  f^  and  f^.  These  equiva¬ 
lent  descr Ipti^.ns  of  a  function  arc  analogous  to  the  study  of  electri¬ 
cal  networks  in  either  the  time  domain  (impulse  and  correlation 
function),  or  frequency  domains  (frequency  response  and  transform 
analyses).  The  appropriate  mathematical  relations  involved  in 
space  filtering  are  described  more  fully  below. 

In  the  general  case,  an  object  has  a  particular  location  in 
space  and  time,  and  also  has  particular  reflection  or  emission 
characteristics  for  electromagnetic  energy.  The  electromagnetic 
energy  is,  for  all  intents  and  purposes,  incoherent  in  the  visible 
and  infrared  spectrums  and  thus  only  the  intensity  is  of  interest. 

The  optical  imaging  system  is  further  assumed  to  be  linear  (dis¬ 
placements  from  center  of  field  and  intensities  are  proportional  in 
object -image) ,  as  is  the  transmissiocA  factor  of  the  nptical  elements: 
i.  e.  ,  they  ire  isoplanatic. 

The  total  optical  processing  from  the  object  plane  to  the  cell 
is  illustrated  in  Figure  40.  A  lens  or  mirror  forms  a  two- 
dimensional  image,  which,  grossly  speaking,  is  linearly  related  to 
the  object.  The  blur  circle  or  aberration  in  the  image  plane  is 
defined  by  the  light  distribution  which  arises  from  a  point  source 
and  is  ajjalogous  to  the  impulse  response  in  electrical  networks. 

The  image  of  an  extended  object,  however,  is  formed  by  the  summa¬ 
tion  of  many  blur  circles  arising  from  the  rnciny  point  sources  in  the 
object.  This  resulting  addition  to  describe  a  function  in  a  linear 


system  is  known  as  the  superposition  integral  and  is  conveniently 
dealt  with  mathematically  by  use  of  the  Fourier  transforms.  There¬ 
fore,  the  object  can  be  expressed  in  terms  of  its  two-dimensional 
space  frequency  content,  the  lens  or  mirror  has  a  space  frequency 
transfer  characteristic  (the  transform  of  the  blur  circle)  or  "pass- 
band”  which  is  an  ability  to  translate  this  object  detail  into  image 
detail,  and  the  resulting  image  has  a  quality  or  space  frequency 
content  which  is  simply  the  product  of  the  object  function  and  lens 
function.  This  dependence  on  the  object-image  function  has  been 
described  by  Duffieux,  Elias,  and  others. 

In  real  space  the  reticle  and  field  stop  combination  multiply 
their  transmission  function  by  that  of  the  image.  The  field  lens 
subsequently  integrates  the  total  light  intensity  which  passes  through 
the  reticle  by  imaging  the  primary  collector.  The  reticle  function 
itself,  however,  must  be  separated  into  its  time  and  space  variables 
if  it  is  simultaneously  chopping  and  scanning.  The  resulting  inte¬ 
grated  space  frequency  multiplication  results  in  a  real  convolution 
which  is  modulated  at  a  time  frequency  depending  on  the  separated 
time  variable  in  the  reticle  function. 

Mathematically,  the  output  of  a  cell  or  detector  c(x,  y,  t) 
is  desired  as  a  function  of:  (1)  the  object  intensity  distribution 
corresponding  to  the  background  or  target,  o(4»  t));  (2)  the  trans¬ 
mission  factcr  of  the  optids,  g(x,  y),  which  corresponds  to  the 
"spread  func.  ion"  or  the  blur  circle  pattern  of  a  point  object  as 
seen  on  the  in  rage  plane;  and  (3)  the  reticle  function,  r(x,  y,  t), 
which  may  be  simultaneously  scanning  and  chopping  the  image. 

The  spread  function,  s(x,  y),  is  the  light  intensity  distribution 
of  a  point  object  in  the  image  plane.  An  extended,  incoherently 
illuminated  object,  o(4»t]),  has  a  resulting  image,  i(x,  y),  composed 
of  the  integral  of  many  point  sources  'n.  the  object  modified  by  the 
transmission  function.  If  unit  magnification  is  assumed,  then 

fr” 

i(x,  y)  =  U  s{x  -  y  -  Tl)  o(t,  ii)  dri  (j) 
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The  light  distributiori  immediately  behind  the  reticle  is  the 
real  space  multiplication  of  the  reticle  and  image.  The  field  lens. 
however5  integrates  this  total  light  intensity  distribution  upon  ^ 
cell  or  detector.  Furthermore,  relative  motion  of  the  reticle  in 
space  imparts  a  .1  space  convolution  to  the  two  space  functions 
when  combined  with  the  integrating  effect  on  the  cell- condenser 
combination.  The  reticle  function  itself  is  separated  into  orthog¬ 
onal  time  and  space  variables  in  a  manner  analogous  to  traveling 
wave  theory.  The  total  light  intensity  falling  on  the  cell  is  the 
product  of  each  individual  differential  area  summed  or  integrated 
over  the  va>rious  reticle  displacements  x\  y',  and  modulated  at 
the  time  frequency  of  the  reticle. 


c(x,  y,  t)  =  r  ^  (t) 


r2(x',  >  )  i{x -X  ,  y -y  )  dx  dy  =  c(t)  c(x.  y)  (2) 


The  functional  relationships  in  Equations  (1)  and  (2)  may  also 
be  related  to  each  other  in  space  frequency  terms  by  employing  the 
two-dimensional  Fourier  transform  of  a  function,  f(x,  y),  defined 


as  follows: 


F(Wx.  Wy)  = 


f(x,  yja’^^^wx^  ^  wy^^dx  dy 


The  resultant  spatial  transform  of  the  cell  output  related  to  all  the 
optical  and  reticle  components  is  then 


wy^  ”  ^2^wx’  wy^  ^^wx’  wy^  ^^wx’  wy^ 


In  summary,  the  total  image  processing  from  the  object  plane 
to  the  detecting  element  can  be  treated  conveniently  with  Fourier 
transforms.  The  light  intensity  modulation  on  the  cell  is  the  prod¬ 
uct  of  the  two-dimensional  transform  or  space  frequency  of  the  object 
the  lens  and  the  reticle  or  field  stop  all  modulated  at  a  time  frequency 
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determined  by  the  separated  time  function  of  the  reticle.  Figure  41 
illustrates  two  examples  of  circular  field  stops  with  rotating  reticles. 

With  stationary  reticles,  the  time  frequency  output  f^  is  equal  to 
the  resulting  transformed  noise  output  f{x,  y),  v  where  v  is 
velocity  in  the  scan  direction. 

In  order  to  estimate  background  interference  in  any  particular 
system,  therefore, the  product  of  the  two  dimensional  noise  spectrum, 
and  the  two  dimensional  reticle  function  are. integrated  at  the  modu¬ 
lating  or  time  frequency  of  the  reticle.  If  for  example,  the  back-  ' 
ground  noise  of  a  3°  field  of  view  tracker,  with  3  milliradian  reticle 
spacinginthe  3.  5  -  5.5p  spectral  band  were  desired,  the  following 
steps  could  be  taken. 

The  reticle  reinforcement  frequency  occurs  at  i/(2  x  3  mr)  or 
160  waves /radian.  The  effective  bandwidth  of  the  aperture  at  this 
frequency  is  approximately  1/Ax  in  the  fx  direction  and  1/Ay  in  the 
fy  direction  or  19  waves /radian  in  each  direction.  The  background 
at  160  waves /radian  is  assumed  while  within  the  space  bandwidths 
of  the  aperture  considered.  The  reticle  has  a  field  of  view 
=  277  X  10*^  steradians  but  it  is  only  50  percent  transmissive. 

The  background  contrast  is  20  percent.  Total  300®K  blackbody 
radiation  is  1.  5  x  10’^  watts /cm^/steradian  of  which  2.  5  percent 
occurs  between  3.  5  and  5.  5p.  The  background  from  Figure  38 
at  l60  waves /radian  for  an  average  distance  to  clouds  of  5  n.  mi. 
is  then 

N(f;^,  f  j  y -  x/l.  5  X  10  ^  X  X  1/2  X  fw/cm^/steradian 

""  ^  160^  V  100  100/  V  j 


lian^ 


I  wave /radian)  =  2  x  10 

The  actual  interfering  noise  irradiance  is 


(w/cm^/steradian)  ^/  (wave /radian)  ‘ 

n  /  \  / 


N(fx.  fy)  R(fx.  fy)  d  fxd  fy 


where  R(fx,  fy)  is  the  Fourier  transform  of  the  reticle  function  or 
approximately  (11 /2)^  (1/Ax)  (1/A  y),  at  the  reinforcing  frequencyi 


r 
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Figure  33.  Radiometric  Sky  Observations 
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Figure  34,  Assumed  Probability  Distribution  of  Cloud  Lengths 
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e  Dimensional  Wiener  Background  Spectra  -  Solar  Reflect: 


Wiener  Background  Spectra  -  Thermal  Emission 


(CONTRAST  IRRADIANCE  LEVEL)  /(WAVE/RADIAN) 


Figure  39»  Predicted  Two  Dimensional  Wiener  Background  Spectra 
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FIELD  LENS 


r,(t)  C(« 

I  I 


Figure  40.  Linear  Optical  Imaging  System 
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PROJECTION  ON  (d 


Transform  of  Parallel  Slit  and  Radical  Reticle 


SUMMARY 


The  investigation  of  infrared  sky  backgrounds  has  been  a 
continuing  project  since  1953.  In  that  year,  at  North  American 
Aviation,  Inc.,  under  Contract  AF  33  (600)-24347  (WADC),  the 
instrument  knov/n  as  a  radiometer-gradient  meter  was  designed, 
fabricated,  and  functionally  tested.  Additionally,  a  high-altitude 
weather  study  was  performed.  Four  quarterly  reports  and  a 
final  report  were  issued,  The  four  quarterly  reports  embody 

the  engineering,  design  and  test  of  the  instrument,  associated 
electronics  and  tlie  low  temperature  black-body  calibrator.  The 
second  quarterly  report  contains  a  section  devoted  to  the  linear 
scanning  motion  of  the  instrument  and  the  difficulties  associated 
with  the  separation  of  the  background  Wiener  spectrum  and  the 
spectra  developed  by  the  distortions  of  a  saw  tooth  or  triangular 
wave  scan.  Other  detaiils  pertinent  to  the  design  of  an  instrument 
to  produce  an  accurate  measure  of  space  frequencies  is  embodied 
in  that  report.  - 

Quarterly  Progress  Report  No.  4  of  AF  33(600)-24347  con¬ 
tains  a  pe  .inent  discussion  on  the  Infrared  Scattering  of  Sky  Light 
based  or.  ere  work  by  R.  Tousey  and  E.  O.  Hulbert.  The  final 
report  fer  .ne  North  American  Aviation  effort  publishes  the  first 
Wiener  spectra  produced  by  the  gradient  meter  in  testing  its 
functional  capability.  '  -alization  of  the  need  for  this  type  of  data 
is  enhanced  by  the  fa».t  tliat  copies  of  this  data,  first  published  in 
1954,  were  requested  of  the  author  as  late  as  1957.  The  final 
report  also  contains  the  results  of  the  high-altitude  weather  research 
which  is  a  very  complete  and  concise  discussion  of  upper  air  meteor¬ 
ological  phenomena.  A  remarkably  complete  and  up-to-date  (1954) 
source  bibliography  on  the  subject  is  included  as  well. 

Infrared  sky  background  investigations  were  continued  by 
Ramo-Wooldridge,  Corp.  ,  beginning  in  1955  under  Contract  AF 
33  (600)-30489.  Three  reports  were  issued  semi-annually,  tracing 
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the  progress  of  the  program.  The  equipment  was  readied  for  field 
use  and  a  company- owned  mobile  laboratory  was  outfitted  for  this 
purpose.  In  addition,  the  radiometer  was  equipped  with  a  set  of 
7  narrow  band  (Turner)  interference  filters  for  making  spectral 
measurements  in  the  region  from  1  to  4.25  microns.  During  this 
program,  the  equipment  was  moved  to  the  Colorado  Springs, 

Colorado  area  for  cooperative  participation  in  the  IRMP/56.  The 
second  semiannual  report  and  the  final  report  for  this  contract 
contain  many  individual  Wiener  spectra  made  at  Mt.  Wilson  at 
6000  feet  altitude  and  also  from  the  roof  of  one  of  the  Ramo- 
Wooldridge  buildings  at  sea  level.  The  final  .  jort  for  this 
contract  contains,  as  an  appendix,  a  reprint  of  a  paper  by  P.  R.  Karr 
of  the  Ramo-Wooldridge  Ccrp.  entitled  ’’Mathematical  Study  of  Back¬ 
ground  ‘Noise’,  "  presented  at  the  second  meeting  of  'die  East  Coast 
IRIS,  26  April  1957.  This  is  a  highly  valuable  paper  because  it 
contains  a  number  of  elements  particularly  concerned  with  circular 
scanning,  which  is  the  method  most  used  in  obtaining  the  gradient 
data.  Some  of  the  equations  are  utilized  in  the  main  body  of  the 
report  to  present  an  example  of  the  space  frequency  transmission 
of  a  simplified  reticle. 

The  background  investigation  work  continued  to  be  performed 

by  Ramo-Wooldridge,  now  a  division  of  Thompson  Ramo  Wooldridge, 

Inc.,  under  Contract  AF  19  (604)-3473,  from  the  Geophysics  Research 

Directorate  of  the  Air  Force  Cambridge  Research  Center.  Four 

reports  have  been  issued  under  this  contract  including  this,  the  final 
(4) 

report.  Report  No.  1  presented  the  first  Wiener  spectra  to  give 
an  indication  of  the  space  frequency  amplitudes  in  the  lead  selenide 
region.  Report  No.  2  is  a  comprehensive  tabulation  of  sky  gradients 
in  the  Key  West  area,  while  No.  3  tabulates  radiometric  measure¬ 
ments  made  at  the  same  time.  The  present  final  report  represents 
the  generalization  of  the  tabulated  data  in  the  two  previous  reports 
and,  in  addition,  presents  a  model  for  background  sky  noise  that 
may  be  of  help  in  predicting  levels  at  the  higher  waves  per  radian 
where  much  interest  is  centered  today. 
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The  equipment  used  since  1953  is  essentially  the  same  as 
designed,  with  a  few  minor  improvements  added  during  vhe  course 
of  time.  It  is  substantial  and  dependable  with  no  instrument 
''headaches".  This  dependability  and  accuracy  of  measurement 
has  enabled  ihe  mass  of  data  to  be  taken.  In  addition,  the  large 
amount  of  gradient  data  taken  at  Key  West  might  have  remained 
unreduced  and  consequently  of  no  use  had  it  not  been  for  the  use 
of  the  computer  facilities  at  Ramo-Wooldridge. 
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